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Introduction 

The  measurement  of  interfacial  tensions  of  fluid  interphases  is  a  problem  of 
strong  current  interest  for  a  broad  range  of  scientific  disciplines  and  tech¬ 
nologies.  As  a  fundamental  parameter,  surface  tension  is  the  single  most 
accessible  experimental  parameter  that  describes  the  thermodynamic  state  and 
structure  of  an  interphase.  The  technology  of  interphases  is  important  to  a 
multitude  of  applications  encompassing  such  diverse  areas  as  coatings,  disper¬ 
sions,  adhesion,  and  membranes  to  name  but  a  few.  It  is  not  surprising  then 
that  a  great  deal  of  effort  has  been  focussed  on  the  development  of  methods  for 
measuring  interfacial  tensions  (see  references  1-3  for  extensive  reviews).  The 
most  versatile  of  the  available  methods  are  those  based  on  the  characterization 
of  axisymmetric  fluid  drop  profiles  (2,4-8). 

The  technique  which  has  found  the  most  widespread  usage  involves  the  deter¬ 
mination  and  analysis  of  pendant  fluid  drop  profiles.  The  shape  of  a  pendant 
fluid  drop  is  specified  in  terms  of  the  cartesian  coordinates  x  and  z  as 
depicted  in  Figure  1,  where  <|)  is  the  angle  between  the  horizontal  axis  and  a 
tangent  to  the  drop  profile  at  position  (x,z).  The  parameter  b  refers  to  the 
radius  of  curvature  at  the  drop  apex  and  s  is  the  contour  length  along  the  pro¬ 
file. 

The  equations  representing  this  profile  can  be  expressed  in  dimensionless 
form  (9)  as 


where  dimensionless  reduced  variables  are  defined  as  X  =  x/c,  Z  =  t\/c,  S  =  s/c, 
and  B  =  b/c. 

The  dimensionless  shape  factor  is  given  by 

B  =  b/c  =  [2] 

where  Ap  is  the  mass  density  difference  between  the  drop  fluid  and  the 
surrounding  medium,  g  is  the  gravitational  constant  (negative  for  pendant  drops) 
and  y  is  the  interfacial  tension. 

Classical  methods  for  the  determination  of  interfacial  tension  rely  on  pho¬ 
tographic  imaging  of  the  drop  profile.  Characteristic  dimensions  of  the  drops 
are  measured  manually  from  enlargements  of  these  photographs,  and  an  estimate  of 
the  interfacial  tension  is  furnished  by  comparison  of  these  dimensions  to  tabu¬ 
lated  theoretical  values  obtained  from  numerical  solution  of  [1]  (1,10). 

Although  simple  in  their  execution,  these  methods  are  limited  in  accuracy  by  the 
necessity  for  specification  of  the  drop  extrema  (i.e.  position  of  apex  and  maxi¬ 
mum  diameter).  In  addition,  the  classical  methods  can  be  quite  tedious,  espe¬ 
cially  for  polymeric  systems  which  require  substantial  equilibration  times. 

The  pendant  drop  techique  has  recently  benefitted  from  a  number  of  advan¬ 
cements  both  in  terms  of  data  acquisition  and  analysis.  The  simplest  of  these 
improvements  makes  use  of  digitizer  palettes  to  efficiently  record  and  store 
profile  coordinates  from  photographic  drop  images  (11-13).  Photographic 
recording  can  be  eliminated  entirely  by  direct  digitization  of  the  drop  image 
with  the  aid  of  either  a  video  frame  grabber  or  a  direct  digital  camera  (14,15). 
There  are  a  number  of  advantages  to  digital  processing  of  drop  images  including 
rapid  acquisition,  facilitated  data  manipulation  and  analysis,  and  the  capabi- 


lity  for  data  storage.  In  addition,  the  techniQues  provide  a  simple  means  to 
detect  the  attainment  of  equilibrium,  a  distinct  advantage  for  viscous  fluids 
such  as  polymer  melts.  This  is  accomplished  by  monitoring  the  time  evolution  of 
the  drop  profile  as  it  asymptotes  towards  its  ^uilibrium  shape. 

The  ability  to  easily  detect  the  entire  drop  profile  has  provided  an  impetus 
for  the  concurrent  development  of  sophisticated  algorithms  for  drop  profile  ana¬ 
lysis  (14-17).  These  methods  either  eliminate  or  minimize  the  requirement  for 
specification  of  extremal  drop  dimensions,  thereby  reducing  the  inherent  sta¬ 
tistical  error.  The  interfacial  tension  can  be  determined  by  performing  a  least 
squares  fit  of  the  segmented  drop  image  to  the  numerical  solution  of  [1].  To 
date  two  methods  for  optimization  of  this  fit  have  appeared;  exhaustive  search 
through  the  shape  parameter  B  (14);  or  a  sophisticated  Newton-Raphson  procedure 
in  conjunction  with  incremental  loading  (16). 

In  this  communication  we  describe  some  further  improvements  in  drop  image 
analysis  including  a  new  shape  comparison  algorithm  for  the  determination  of 
interfacial  tensions  from  axisymmetric  fluid  drop  profiles.  The  algorithm 
employs  robust  shape  comparison  routines  (18-20)  which  are  resistant  to  outlying 
points  that  may  result  from  inaccuracies  in  the  computerized  drop  profile 
discrimination  procedures.  Furthermore,  the  algorithm  is  designed  to  be  imple¬ 
mented  on  a  microcomputer,  thus  enabling  routine  automation  of  absolute  boundary 
tension  measurements. 


Experimental 

A  block  diagram  of  the  experimental  apparatus  for  boundary  tension  measure¬ 
ments  appears  in  Figure  2.  Fluid  drops  are  formed  at  the  tip  of  either  a 
stainless  steel  syringe  or  glass  capillary  tube  which  is  inserted  through  a  sep- 


turn  into  a  quartz  cuvette.  The  cuvette  is  mounted  within  a  cell  with  provisions 
for  both  atmosphere  and  temperature  control  (±1*C).  The  optical  system  consists 
of  a  Unitron  microscope  (zoom  stereo  trinocular-ZST)  coupled  to  an  RCA  model  TC 
1005  video  camera.  The  optics  are  focussed  by  optimizing  the  video  image  of  a 
reticle  containing  a  finely  ruled  grid  that  is  placed  at  the  drop  location".  In 
addition  to  furnishing  a  well  defined  means  for  focussing,  the  grid  facilitates 
the  detection  of  any  image  barreling  or  other  distortions  that  may  arise  from 
the  optics,  and  provides  a  direct  calibration  of  both  the  vertical  and  horizon¬ 
tal  magnification  factors  inherent  to  the  optics  and  camera  systems. 

The  signal  from  the  video  camera  is  fed  to  a  Tecmar  Video  Van  Gogh  board, 
which  performs  the  frame  grabbing  and  image  digitization.  In  the  high  resolu¬ 
tion  mode,  the  image  is  comprised  of  a  251  by  255  pixel  array  with  256  gray 
levels.  The  image  analysis  software  routines  are  coded  in  machine  language  and 
Fortran  and  are  executed  on  an  IBM  PC. 

Drop  Profile  Analysis 
A.  Segmentation 

The  process  of  drop  profile  discrimination  (i.e.  segmentation)  is  carried 
out  by  applying  global  thresholding  to  the  drop  image.  Global  thresholding 
reduces  the  256  gray  level  image  to  a  binary  image  by  judicious  choice  of  an 
appropriate  pixel  intensity  magnitude  or  threshold.  In  order  to  discriminate 
the  drop  interface,  the  proper  threshold  magnitude  must  correspond  to  an  inten¬ 
sity  level  intermediate  between  that  of  the  drop  and  its  surrounding  medium. 

The  optimal  threshold  value  is  chosen  with  the  aid  of  a  grey  level  histogram. 
Castleman  (21)  has  presented  arguments  suggesting  that  the  optimal  threshold 
value  is  located  at  the  intermediate  minimum  in  the  bimodal  gray  level 


histogram.  Alternative  segmentation  techniques,  such  as  adaptive  thresholding 
(21)  or  gradient  magnitude  techniques  (21-26)  could  also  be  adopted.  ^  The 
algebraic  operations  for  these  alternative  methods  require  considerable  execu¬ 
tion  time  however  without  substantial  improvement  in  accuracy  (under  proper 
lighting  conditions)  as  suggested  by  PavTidis  (22).  Figure  3  shows  typical 
digital  images  of  a  pendant  drop  before  thresholding  (left),  after  thresholding 
(center),  and  the  resultant  segmented  drop  profile  (right).  A  simple  test  for 
successful  edge  detection  is  superposition  of  this  profile  on  the  original  drop 
image . 


B.  Piece-Wise  Rotational ly  Resistant  Smoothing 

The  resolution  of  the  segmented  drop  profile  is  limited  by  the  resolution  of 
the  frame  grabber  as  can  be  seen  from  Fig.  3  (right).  The  situation  can  be 
improved  by  the  use  of  higher  resolution  frame  grabbers,  but  this  leads  to 
increases  in  both  the  cost  and  execution  time  of  the  instrument.  An  alternative 
approach  to  increase  resolution  is  to  convert  the  integer  pixel  values  of  the 
digitized  drop  profile  coordinates  into  the  appropriate  floating  point  decimal 
numbers.  This  can  be  accomolished  by  the  application  of  local  polynomial 
regression  smoothing  methods.  The  smoothing  operations  are  applied  piece-wise 
(i.e.  point-by-point  replacement)  along  the  entire  drop  profile  contour  to 
complete  the  smoothing  procedure.  Smoothing  parameter  options  include  the  order 
of  the  polynomial  smoothing  (e.g.  quadratic  or  cubic),  the  number  of  contiguous 
data  points  in  the  local  target  group,  and  the  number  of  overall  smoothing 
cycles.  Smoothing  techniques  are  often  subject  to  overshoot  in  regions  of  steep 
slope  such  as  are  encountered  in  the  near-circular  shapes  of  pendant  fluid 
drops.  In  order  to  avoid  overshoot,  a  rotationally  resistant  smoothing  routine 


was  developed. 


The  smoothing  algorithm  begins  by  calculating  the  local  slope  of  an  odd 


number  of  contiguous  points  (e.g.  21)  by  a  principle  components  fit  (i.e.  linear 
regression).  These  "target"  points  are  then  rotated  using  the  resultant  coor¬ 
dinate  transformation  matrix  such  that  the  slope  is  zero.  In  addition,  the 
points  may  be  translated  to  the  origin  in  order  to  normalize  the  magnitudes  of 
the  position  coordinates.  After  rotation  and  translation,  smoothing  is 
accomplished  by  performing  a  least  squares  or  robust  polynomial  fit  in  the  new 
coordinate  system.  The  smoothed  profile  coordinate  (i.e.  midpoint  of  each 
target  group)  in  the  original  coordinate  system  is  then  obtained  by  a  second 
application  of  the  coordinate  transformation  matrix  that  returns  the  profile 
segment  to  its  original  position. 

C .  Shape  Fitting 

The  process  of  fitting  an  experimental  drop  profile  to  the  theoretical  shape 
predicted  by  [1]  involves  a  five  parameter  optimization.  Three  parameters  are 
required  for  alignment  of  the  drop.  They  are  a  rotation  angle,  an  x  translation 
and  a  z  translation.  Drop  rotation  is  required  in  order  to  align  the  camera 
system  to  the  vertical  axis  of  symmetry  of  the  drop.  The  final  two  parameters 
required  are  the  scale  or  magnification  factor  of  the  drop  (/c;  see  [1]),  and 
the  shape  factor,  B.  The  axis  of  symmetry  is  determined  by  minimizing  the  sum 
of  squared  differences  between  the  distances  from  corresponding  "symmetric" 
points  on  each  side  of  the  drop  to  the  axis  of  symmetry.  In  practice,  this 
algorithm  is  found  to  be  insensitive  to  drop  symmetry  due  to  the  near-spherical 
shape  of  the  drop.  It  is  preferable  to  eliminate  the  need  for  rotation  alto¬ 
gether  by  prealigning  the  vertical  camera  axis  with  the  direction  of  gravity. 

The  comparison  of  shapes  by  robust  methods  requires  establishment  of  a 
reference  point  so  that  specific  points  of  comparison  may  be  selected.  There 


are  two  natural  choices  for  this  reference  point:  the  drop  apex,  or  the  drop 
"center"  defined  as  the  intersection  of  the  vertical  axis  of  symmetry  and  the 
horizontal  maximal  drop  diameter.  Points  for  comparison  from  the  experimental 
profile  are  chosen  at  predetermined  intervals  along  the  contour  of  the  drop  pro¬ 
file.  Corresponding  points  from  the  theoretical  profile  are  selected  by 
requiring  that  the  line  segment,  connecting  the  point  of  interest  and  reference 
point  on  the  theoretical  profile,  have  the  same  slope  as  the  corresponding  line 
segment  on  the  experimental  profile.  Points  on  both  sides  of  the  drop  are 
employed  in  order  to  ensure  symmetry. 

A  first  estimate  of  the  dimensionless  shape  parameter,  Bq,  is  obtained  from 
the  approximate  relation  (17) 

Bq  =  exp[-6. 70905  +  15.3002W  -  16.44709w2  +  9,92425w3  -  2.58503W^] 

[2] 

where  W  is  the  ratio  of  the  local  droo  c^ameter  (for  a  plane  located  at  a 
distance  from  the  drop  apex  which  is  equal  in  magnitude  to  the  maximum  droo 
diameter)  to  the  maximum  drop  diameter. 

The  optimal  B  value  is  obtained  by  exhaustive  search  in  the  vicinity  of  Bq- 
This  optimization  is  accomplished  by  robust  shape  comoarison  routines  developed 
by  Siegel  and  co-workers  (18-20).  As  in  all  regression  problems,  the  process 
involves  fitting  N  points  (x,j,yi)  from  an  experimental  profile  to  N  homologous 
points  (u,j,v^)  on  a  theoretical  curve.  The  theoretical  points  must  be  rotated 
by  an  angle  8,  translated  by  a  vector  (a, 8)  and  scaled  by  a  factor  t  in  order  to 
affect  this  comparison.  The  transformed  theoretical  coordinates 


are  then  compared  to  (x^.y^)  for  each  value  of  the  shape  factor,  B.  The  value 
of  the  shape  factor  which  yields  the  minimum  overall  error  provides  the  optimal 
fit.  The  interfacial  tension  is  then  obtained  from  the  associated  optimal 
scaling  factor  t,  by  recognizing  that  t  =  i^c  =  (Apg/y)  1/2, 

In  least  squares  regression  methods,  the  values  of  all  of  the  shape  parame¬ 
ters  (i.e.,  T,e,o,/3)  must  be  chosen  simultaneously  in  order  to  minimize  the  sum 
of  squared  residuals 

N  '2  '2 

sum  =  I  [(X,.  -  u.)  +  {y,.  -  v^)  ]  [4] 

i  =  1 

In  contrast,  with  robust  shape  comparison  methods,  each  of  the  optimal  paramete 
values  may  be  evaluated  independently.  In  the  case  of  the  rotation  and  magnifi 
cation  variables,  this  is  accomplished  using  the  concept  of  repeated  medians 
(18-20)  as  represented  by  the  relations. 

T  =  med  {med  {t.  .) )  [5] 

i  j  ^ 


where 


9 


where  is  the  rotation  required  for  the  vector  from  point  i  to  point  j  in  the 
theoretical  shape  to  have  the  same  attitude  as  the  homologous  vector  ^in  the 
experimental  profile. 

The  translation  parameters  are  calculated. from  the  simple  (non-repeated) 
medians : 


a  =  med  {x.  -  T[u.cose  -  v.sin©])  [8] 

i 

3  =  med  {y.  -  T[u.sine  +  v.cosS]}  [9] 

i 

The  advantages  of  double  median  robust  techniques  over  traditional  least 
square  regression  methods  are  discussed  in  references  (18-20). 

For  each  B  value  assumed,  the  optimal  values  of  t,  B,  a,  and  are  specified 
directly  by  these  relations.  The  overall  optimization  routine  is  therefore 
reduced  to  an  exhaustive  search  on  the  shape  parameter  B.  The  ultimate  best 
fit,  and  associated  interfacial  tension,  is  determined  from  the  B  value  which 
produces  the  minimum  overall  error  or  deviation  from  the  experimental  profile. 
Three  distinct  types  of  error  may  be  calculated  for  this  purpose. 
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The  term  FACTOR  normalizes  the  total  error  to  the  overall  size  of  the 
object.  TENSE  is  a  measure  of  the  "stress"  or  "tension"  required  in  "deforming" 
the  theoretical  profile  so  as  to  fit  the  experiment,  and  in  some  cases  is  a  more 
effective  measure  of  shape  differences  than  is  the  simple  root  mean  square 
error,  SR.  While  the  root  mean  square  statistic  measures  the  magnitude  of  the 
differences,  the  tension  measures  how  quickly  this  difference  is  changing.  The 
error  SRR  is  defined  in  a  fashion  to  be  more  resistant  to  outlying  points  than 
the  standard  sum  of  squared  residual,  SR. 


Oj.  Results  of  Simulations 

The  performance  of  the  shape  comparison  algorithm  was  examined  by  applying 
it  to  simulated  pendant  drop  profiles  generated  by  a  Runge-Kutta  numerical  solu¬ 
tion  of  [1].  A  typical  theoretical  profile  is  shown  as  the  solid  line  in  Figure 
4,  for  a  dimensionless  shape  parameter  of  B  =  0.6.  In  order  to  simulate  the 
actual  conditions  encountered  with  the  pendant  drop  apparatus,  several  sources 
of  error  were  introduced.  The  figure  at  left  shows  the  same  profile  after 
introduction  of  Gaussian  random  noise  with  an  average  magnitude  of  error 
corresponding  to  one  pixel  element  (points).  The  figure  at  right  shows  the  ori¬ 
ginal  theoretical  profile  (solid  line)  along  with  a  profile  (points)  that  has 
been  rotated  by  one  degree  to  simulate  misalignment  of  the  camera  and  drop  pro¬ 
file  coordinate  systems.  This  latter  profile  was  generated  using  the  same 
Gaussian  error  distribution  as  for  the  former  profile,  but  also  includes  addi- 


tional  non-random  errors  that  were  added  arbitrarily.  Non-random  errors  were 
added  in  order  to  simulate  the  possible  effects  of  improper  edge  detection  that 
might  occur  for  images  of  low  contrast. 

The  smoothing  algorithm  that  incorporates  translation  of  the  target  points 
to  the  origin  as  well  as  rotation  provided  the  best  performance,  especially  in 
regions  near  the  maximum  drop  diameter.  This  result  may  be  attributed  to 
the  decreased  magnitudes  of  the  drop  profile  coordinates  after  they  are  shifted 
to  the  origin.  A  drawback  of  translation  to  the  origin  is  the  increase  in 
computation  time. 

The  qualitative  performance  of  the  smoothing  algorithms  is  demonstrated  in 
Figure  5.  It  is  apparent  that  the  presence  of  non-random  error  (righthand 
figure)  induces  considerable  distortion  into  the  smoothed  drop  profile  in  the 
vicinity  of  the  outlying  point.  In  principle,  these  effects  could  be  reduced  by 
the  application  of  robust  rather  than  least  squares  smoothing  routines. 

The  quantitative  performance  of  the  shape  analysis  and  smoothing  algorithms 
was  assessed  by  applying  them  to  the  simulated  data.  When  the  profile  was 
smoothed  prior  to  shape  analysis,  comparable  results  were  obtained  for  both 
choices  of  reference  point  (i.e.  apex  or  center).  For  unsmoothed  profiles  sub¬ 
ject  to  considerable  noise,  however,  use  of  the  drop  center  yielded  superior 
results.  This  latter  observation  may  be  easily  explained  by  the  fact  that  the 
slopes  of  line  segments  that  define  the  comparison  points  are  less  sensitive  to 
error  for  the  center  reference  point  than  for  the  apex  reference  point.  For 
highly  elongated  drops,  however,  the  apex  becomes  the  reference  point  of  choice, 
since  definition  of  the  position  of  maximum  diameter  is  problematic.  In  fact, 
the  apex  must  be  chosen  as  the  reference  point  for  drops  that  do  not  exhibit  a 


maximum  diameter. 
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The  improvement  in  performance  that  results  from  smoothing  is  illustrated  in 
Figure  6,  where  the  error  SR  is  plotted  as  a  function  of  the  dimensionless  shape 
parameter.  The  leftmost  figure  presents  results  of  the  analysis  applied  to  the 
simulated  drop  profile  (B  =  0.6)  with  random  Gaussian  error.  Without 
smoothing,  the  minimum  error  corresponds  to  a  shape  parameter  of  Bgxp  =  0.597. 
The  deviation  from  the  simulated  profile  decreases  (B^^p  =  0.598)  after  applica¬ 
tion  of  rotationally  resistant  smoothing.  Similar  results  are  obtained  from 
analysis  of  the  profile  with  non-random  errors  and  rotation  (Fig.  6,  right).  In 
this  case  the  magnitudes  of  deviation  from  the  input  value  (B  =  0.6)  are  larger 
as  might  be  expected.  Analysis  of  simulated  drop  profiles  for  B  =  0.5  yielded 
best  fits  of  Bgxp  =  0.494  (without  smoothing),  Bgxp  =  0.501  (with  rotationally 
resistant  smoothing),  and  Bgxp  =  0.500  (rotationally  resistant  smoothing  which 
incorporates  translation  to  the  origin),  again  illustrating  an  improvement  in 
precision  upon  smoothing. 

The  actual  precision  of  the  shape  analysis  algorithm  will  depend  on  the 
shjpe  parameter  B,  and  the  number  of  data  points  selected  for  comparison.  The 
magnitude  of  the  error  can  be  reduced  by  smoothing  and  is  a  function  of  the 
order  chosen  for  the  polynomial  smoothing  function,  the  number  of  target  points 
for  smoothing,  and  the  total  number  of  smoothing  cycles.  Each  improvement  in 
precision  however  is  necessarily  accompanied  by  an  associated  increase  in  com¬ 
putation  time. 

In  most  simulations,  the  error  curves  (e.g.  Fig.  6)  exhibited  discon¬ 
tinuities  near  the  optimal  values  of  B.  In  a  number  of  these  simulations,  the 
least  squares  algorithm  of  Rotenberg  et  al.  (11,16)  failed  to  converge  when 
applied  to  the  data.  The  occurrence  of  this  discontinuity  in  error  may  also  be 
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the  underlying  reason  that  their  algorithm  often  fails  to  converge  if  a  non-zero 
initial  shape  parameter  estimate  (i.e.  Bq)  is  used,  since  the  algorithm  requires 
a  continuous  error  function  for  the  Newton-Raphson  method  to  function  properly. 

Experimental  Application  to  the  Surface  Tension  of  Glycerin 

The  shape  analysis  algorithm  was  applied  to  the  determination  of  the 
surface  tension  of  glycerin  in  order  to  assess  its  performance  under  normal 
experimental  conditions.  The  density  difference  Ap  for  glycerin  was  estimated 
from  density  and  vapor  pressure  data  tabulated  in  Perry's  handbook  (27).  For 
glycerin  at  room  temperature  (24.2  ±0.2*C)  the  vapor/liquid  density  difference 
was  calculated  to  be  Ap  =  1.1745  g/cm3.  The  results  obtained  from  analysis  of 
four  pendant  drop  experiments  gave  a  value  of  62.6  ±0.3  dynes/cm  for  the  room 
temperature  surface  tension  of  glycerin.  This  compares  well  with  the  literature 
value  of  63.4  dynes/cm  (28)  reported  for  a  slightly  lower  temperature  (20*C). 

The  analysis  of  these  drops  was  performed  using  21  data  points  for  the  shape 
comparison.  The  correspondence  of  the  fitted  profile  to  these  data  points  and 
to  the  original  digitized  drop  profile  is  excellent  as  shown  in  Fig.  7. 

The  precision  of  the  experiment,  ca.  0.5%,  is  quite  remarkable  considering 
that  only  21  data  points  were  used  and  that  the  comparison  algorithms  described 
herein  were  designed  specifically  for  execution  on  a  microcomputer.  More 
sophisticated  algorithms  based  upon  extensions  of  these  same  concepts  of  shape 
analysis  can  easily  be  developed  for  implementation  with  more  powerful  main 
frame  computers.  Such  improvements  will  undoubtedly  increase  the  precision  of 
future  interfacial  tension  measurements.  In  addition,  improved  computation 
rates  may  enable  the  future  study  of  transient  interfacial  phenomena. 

f 


Summary 

A  novel  method  for  the  analysis  of  interfacial  tension  from  the  shape  of 
pendant  fluid  drops  is  presented.  The  experimental  apparatus  employs  digital 
video  image  processing  techniques  to  record  drop  images  and  to  segment  the  fluid 
drop  profile.  An  associated  drop  profile  shape  analysis  algorithm  is  developed 
based  upon  robust  shape  comparison  methods  that  utilize  repeated  median  con¬ 
cepts.  The  complete  algorithm  is  designed  specifically  to  be  executable  on  a 
microcomputer.  The  performance  of  the  algorithm  is  illustrated  for  simulated 
pendant  fluid  drop  profiles  and  pendant  drops  of  glycerin.  In  the  latter 
case,  the  experimental  precision  in  the  resultant  surface  tension  is  found  to  be 
of  the  order  of  0.5%. 
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Figure  Captions 


Fig.  1  Experimental  geometry  of  pendant  drop  experiment.  * 

Fig.  2  Block  diagram  of  pendant  drop  apparatus. 

Fig.  3  Typical  digitized  drop  images 

Original  grey  level  image  (left); 

Thresholded  binary  image  (center); 

Segmented  drop  profile  (right). 

Fig.  4  Simulated  pendant  drop  profile.  (Dimensionless  shape  parameter  B  * 


Runge-Kutta  solution  of  [1]  (solid  lines);  points  are  simulated 
integer  number  profiles  with  Gaussian  random  error  (left);  with 
Gaussian  random  error,  non-random  error,  and  a  one  degree  counter¬ 
clockwise  rotation  (right). 

Fig.  5  Effect  of  rotational ly  resistant  smoothing  with  translation  to  the 
origin . 

Points  are  profiles  with  Gaussian  random  error  (left);  with  Gaussian 
random  error,  non-random  error,  and  rotation  (right);  Solid  lines  are 
the  results  of  smoothing. 

Fig.  6  Sum  of  squared  residuals  as  a  function  of  the  dimensionless  shape  para¬ 
meter  for  profiles  simulated  for  B  =  0.6.  Circles:  analysis  without 
smoothing;  Open  squares:  analysis  with  rotationally  resistant 
smoothing;  Filled  squares:  analysis  with  rotationally  resistant 
smoothing  incorporating  translation  to  the  origin. 

Analyis  of  profile  with  Gaussian  random  error  (left); 

Analysis  of  profile  with  Gaussian  random  error,  non-random  error  and 


rotation  (right). 

Digitized  pendant  drop  images  of  glycerin  at  room  temperature  (24.2®C) 
Solid  lines  are  the  theoretical  profiles  obtained  by  application 


of  the  algorithm. 

Data  points  are  the  selected  experimental  data  points  (left); 
Data  points  are  the  original  segmented  experimental  drop  profile 
(right) . 
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